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The Praesepe cluster contains a number of S Set and 7 Dor pulsators. Asteroseismology of cluster stars is simplified by 
the common distance, age and stellar abundances. Since asteroseismology requires a large number of known frequencies, 
the small pulsation amplitudes of these stars require space satellite campaigns. The present study utilizes photometric 
MOST satellite measurements in order to determine the pulsation frequencies of two evolved (EP Cnc, BT Cnc) and two 
main-sequence (BS Cnc, HD 73872) 5 Set stars in the Praesepe cluster. The frequency analysis of the 2008 and 2009 data 
detected up to 34 frequencies per star with most amplitudes in the submillimag range. In BS Cnc, two modes showed 
strong amplitude variability between 2008 and 2009. The frequencies ranged from 0.76 to 41.7 cd _1 . After considering 
the different evolutionary states and mean stellar densities of these four stars, the differences and large ranges in frequency 
remain. 
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1 Introduction 

The majority of stars of spectral type A and F pulsate 
with pressure (p) and/or gravity (g) modes. The stars on 
and near the main sequence with dominant pressure modes 
are known as S Set stars, while the (usually cooler) stars 
with dominant gravity modes are known as 7 Dor stars. 
Typical periods of S Set stars range from 20 minutes to 8 
hours for the dominant modes. Recent campaigns have re- 
vealed hundreds of simultaneously excited pulsation modes 
in some of these stars. The vast majority of the pulsations 
are low-amplitude nonradial modes. The new space- satellite 
Kepler measurements have clearly confirmed the previous 
discovery (e.g., Handler 2009) that many of the A/F star pul- 
sators near the main sequence are hybrid pulsators showing 
both pressure and longer-period gravity modes. 

Besides the 'standard' low-amplitude multimode pul- 
sators, which pulsate with hundreds of nonradial pulsation 
modes, a S Set subgroup is usually recognized: the HADS 

* Based on data from the MOST satellite, a Canadian Space Agency 
mission, jointly operated by Dynacon Inc., the University of Toronto Insti- 
tute for Aerospace Studies and the University of British Columbia with the 
assistance of the University of Vienna. 
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(high-amplitude Delta Scuti) stars are slow rotators with one 
or more dominant radial modes. The commonly used defi- 
nition of HADS stars in terms of peak-to-peak amplitudes 
in excess of 0.3 mag is not meaningful in the case of mul- 
timode pulsators. Also, recent measurements have shown 
that slowly rotating stars with one or two dominant radial 
modes, which are the characteristic of HADS stars, may 
even have smaller amplitudes than this artificial limit, e.g., 
the Kepler star KIC 9700322 (Breger et al. 2011). The 
sizes of the amplitudes and the determination of whether the 
dominant modes are radial or nonradial are a strong function 
of the rotational velocity without an absolute limiting value 
(Breger 2000). Consequently, it is difficult to define (rather 
than just describe) the HADS subgroup. 

Several variability surveys in the Praesepe cluster have 
been undertaken (e.g., Breger 1970; Frandsen et al. 2001, 
Dall et al. 2002) and have revealed a number of low-ampli- 
tude 5 Set stars. The pulsations of these stars are of inter- 
est, because as cluster members they share similar prop- 
erties such as chemical abundance and their distances are 
accurately known. This simplifies the asteroseismic model- 
ing. In this paper we examine the pulsation properties of 
four pulsators studied with the MOST satellite. Here we note 
that MOST satellite measurements were already successfully 
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used in the study of the spectroscopic binary HD 73709 by 
Pribulla et al. (2009, 2010). 

In this paper we present the new results for four pulsat- 
ing members of the Praesepe cluster. More detailed refer- 
ences to previous investigations of these stars are given at 
the end of the discussion of each star together with a com- 
parison of the results. 

2 MOST observations and reductions 

MOST (Micro variability and Oscillations of STars, Walker 
et al. 2003) is a small Canadian spacecraft with a 15 -cm 
telescope designed for photometry of bright stars. Due to 
the polar dusk-dawn orbit, continuous photometry for about 
one month can be obtained for stars close to the ecliptic. 
The nearby open cluster Praesepe was observed from 2008 
January 24 to February 25, and from 2009 February 12 to 
March 1 . MOST employs a custom broadband filter with a 
wavelength range 350 to 750 nm. 

The Direct Imaging (DI) and guide star data were used 
for this investigation. The data were reduced using the data 
reduction developed by Huber & Reegen (2008) and by 
Hareter (2008). The basic function of these data reduction 
programs is based on the decorrelation technique for MOST 
main targets (Reegen et al. 2006), modified to meet the re- 
quirements of the particular data formats. The DI data are 
available as 20x20 pixel images centered on the particular 
star. The guide star data are a compressed version of this 
format, i.e., the intensities of the pixels are added, where 
the mean background is subtracted on-board. Each mea- 
surement consists of this sum and the corresponding inte- 
gration time only. (In contrast to the DI data, no images of 
the guide stars are available). In spite of this on-board back- 
ground subtraction, the data still contain heavy stray light 
which requires additional data reduction. 

The customary three- star photometry reduction scheme 
cannot be applied, since the stray light is distributed very 
inhomogeneously on the CCD chip. The mean background 
value, which is subtracted on-board, is not stored due to 
downlink limitations; hence the true raw data cannot be 
reconstructed. To improve the accuracy of the data, it is 
possible to use the decorrelation technique, as done by the 
previously mentioned data reduction techniques. Here con- 
stant comparison stars are used instead of background pix- 
els, whose individual values are not available. 

During 2008, four S Set stars (BT Cnc, EP Cnc, BS Cnc 
and HD 73872) could be observed. The measurements, con- 
sisting of about 34250 individual measurements for each 
star, covered 32d. During 2009, EP Cnc and BS Cnc were 
measured again with about 17000 measurements per star, 
covering 17d. 

3 Frequency analyses 

The frequency analysis was conducted using the 
PERIOD0 4 software package (Lenz & Breger 2005) 



which utilizes single-frequency Fourier as well as multi- 
frequency least- squares algorithms. The latter technique 
fits a number of simultaneous sinusoidal variations in the 
magnitude domain and does not rely on prewhitening. The 
multifrequency solution obtained in this manner needs to 
be examined in detail in order to (i) determine the statistical 
significance of each detected frequency, and (ii) separate 
the frequencies intrinsic to the star from the instrumental 
effects caused by the MOST satellite. 

The statistical significance of the detected frequencies 
was determined from the amplitude signal/noise (S/N) ra- 
tio, which was calculated in the following manner: for each 
possible new frequency, /(test), the amplitude of /(test) 
was derived from a multifrequency solution with /(test) 
included as one of the frequencies. The noise was deter- 
mined from the best multifrequency solution using all the 
detected significant frequencies excluding /(test). This 
solution was then prewhitened and the average noise calcu- 
lated from the amplitude spectrum using a 5 cd _1 range cen- 
tered around the frequency under consideration. This made 
it possible to derive an amplitude signal/noise ratio with- 
out biassing the outcome by assuming the significance of 
/(test). Furthermore, where data were available from both 
2008 and 2009, the signal/noise ratio was determined from a 
common solution. This guards against any hypothetical spu- 
rious peaks occurring in only a single data set. Following 
Breger et al. (1993) and Kuschnig et al. (1997), a frequency 
was regarded as significant if the amplitude signal/noise ra- 
tio was 4.0 or greater. The final residuals between the data 
and the fits were between 2.3 and 3.1 millimag per single 
measurement. 

We regard our chosen method to include or reject fre- 
quency peaks, based on their amplitude signal/noise ratio, 
as conservative. The chief remaining uncertainty concerns 
any unrecognized systematic instrumental (orbital) effects. 

The remaining stray-light artifacts are far from sinu- 
soidal and not of constant amplitude in each cycle. Hence, 
the stray light frequencies cannot simply be prewhitened. 
A further difficulty arises from the outlier rejection, which 
are most abundant at the phases where the stray light is 
strongest. Moreover, there is also a modulation of the stray 
light with one cycle per day. This leads to a rather compli- 
cated aliasing pattern (spectral window), where the alias- 
ing frequencies are f or b, 2f or b, surrounded by aliasing with 
the modulation frequencies, which is shown in Fig. [2] top 
panel. These pattern also appears in the Fourier transform 
of a MOST data set, not only in the spectral window. Fur- 
ther known artifacts of the MOST photometry can be found 
at 1, 2, 3 cd _1 , usually with decreasing amplitudes. 

4 BT Cnc (= 38 Cnc = HD 73575 = KW 204) 

The data of the star BT Cnc are typical of all four stars. 
Consequently, we can use the 2008 results for BT Cnc to 
illustrate our analysis and techniques. Fig. [T] shows a small 
part of the light curve together with the fit derived below 
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Fig. 1 Sample light curve of BT Cnc covering two-thirds of a day. The continuous curve is the fit derived in this paper. 



and is typical of BT Cnc and the other stars. The small gaps 
in the coverage and the small systematic deviations between 
the fit and the observations are caused by instrumental prob- 
lems of the satellite. 

Table 1 Pulsation frequencies of BT Cnc for 2008 



Frequency Amplitude S/N 

cd _1 millimag ratio 



31 
24 
19 
9.8 
9.5 
9.5 
7.2 
8.0 
7.7 
7.1 
5.5 
6.5 
5.4 
5.5 
5.4 
4.8 
4.8 
4.0 
9.6 
6.4 
6.5 
6.5 



As described in the previous section, the regular fre- 
quencies present in the data were analyzed with the 



PERIOD 4 package. Because of the presence of both stel- 
lar and orbital peaks, we illustrate the intermediate and final 
results in Fig.[2| where panels (a) to (f) are discussed below: 

(a) The top panel shows the spectral window with the 
value of amplitude 1 at cd _1 off scale to illustrate the 
aliasing occurring near the orbital frequency of 14.2 cd _1 
and its multiples. 

(b) This shows the Fourier amplitude spectrum from to 
50 cd _1 . Three frequencies marked as "fl", "f2", and "f3", 
are dominant. 

(c) If we inspect the low-amplitude region (with fi to f3 
off- scale), we see a large number of additional frequencies 
in the 0-50 cd _1 range. 

(d) The number of frequency peaks is severely reduced 
once we prewhiten fi to f3. In fact, except for orbital fre- 
quencies, which are multiples of 14.19 cd _1 , no significant 
peaks remain at frequencies higher than 20 cd _1 . 

(e) However, this still leaves a number of additional sig- 
nificant peaks in the 0-20 cd _1 range. These peaks are 
stellar as well as orbital frequencies, artifacts and aliases. 

(f) 22 statistically significant frequencies are judged by 
us to originate in the star. After we prewhiten these 22 fre- 
quencies, several low-amplitude peaks remain which we 
suspect to be mainly nonstellar because they are close to 
the frequencies of known artifacts. The next stellar peak is 
found at 7.71 cd _1 , but its amplitude signal/noise of 3.7 is 
below our significance limit. Consequently, the search for 
stellar frequencies is stopped at this point. 

The 22 significant stellar frequencies are listed in Table 
[T] The average measurement fits the solution to ± 2.5 mil- 
limag. This also holds for the section of the data shown in 
Fig. [T] The listed frequencies should be free of orbital ef- 
fects with one exception: f 4 at 8.166 cd _1 is separated by 
only 1/T from a suspected orbital alias frequency at 8.199 
cd _1 , where T is the length of observations. A misidenti- 
fication should not be a problem because of the relatively 
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Fig. 2 Amplitude spectra of BT Cnc. The spectral win- 
dow is presented at the top. Panel (c) shows the low- 
amplitude peaks of the 2008 data (with "fl" to "f3" off 
scale). Panel (d) shows that most peaks in panel (c) disap- 
pear after prewhitening the three dominant frequencies. At 
frequencies above 20 cd _1 , only the orbital frequencies of 
the satellite are left. Panels (e) and (f) demonstrate that after 
prewhitening the 22 detected stellar frequencies, most of the 
remaining peaks are instrumental artifacts (marked by 'a'). 



ond and third frequencies reported by Breger (1980), which 
both are 1 cd _1 aliases of frequencies found in the present 
paper. 

5 EP Cnc (= HD 73819 = KW 348) 

For EP Cnc we have detected 20 stellar frequencies, which 
are listed in Table [2] 



Table 2 

2009 



Pulsation frequencies of EP Cnc for 2008 and 
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Amplitude 
millimag 


S/N 
ratio 
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^ The formal errors for the listed frequencies are as small as 

0.0001 cd _1 . However, since annual alias values cannot be 

excluded, we list the value of the annual alias. 

2 ) The annual alias value of fis at 11.5043 cd _1 fits equally well. 

The alternate value can be identified as a combination frequency 

(f 4 +fn). 



large amplitude of f±. Furthermore, after prewhitening 
we are left with a small peak at 8.199 cd _1 , the expected 
orbital alias frequency. 

How do the new frequencies compare with those found 
in previous, smaller studies (Breger 1973; Guerrero, Man- 
tegazza & Scardia 1997; Breger 1980; Kim & Lee 1995; 
Freyhammer, Larsen & Petersen 1997; Hernandez et al. 
1998; Pena et al. 1998)? All studies have detected the domi- 
nant frequency of 9.78 cd _1 , while the 7.02 cd _1 frequency 
usually was also seen. We interpret the differences in the 
other reported values for suspected frequencies to be due to 
the relatively small data sets and 1 cd _1 aliasing of earth- 
based telescope measurements. Examples may be the sec- 



We also calculated individual amplitudes for 2008 and 
2009 separately. The variation of amplitudes between 2008 
and 2009 was small with an average deviation of 0.10 mil- 
limag, except for the low frequency where the amplitude 
changed from 0.9 to 1.2 millimag. 

The results are in excellent agreement with earlier mea- 
surements. The variability of EP Cnc had been discovered 
by Martin & Hube (1989) with an estimate of a 5.46 cd" 1 
frequency. The multisite campaign by the Delta Scuti Net- 
work determined three frequencies (Breger et al. 1994), 
which were identical to the three modes with the highest 
amplitudes in the present data sets. 
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6 BS Cnc (HD 73450 = KW 154) 

For BS Cnc we have detected 20 stellar frequencies, which 
are listed in Table [3] Due to the strong amplitude variabil- 
ity associated with some of the detected frequencies, we list 
the amplitudes from 2008 and 2009 separately. However, 
the detection significance, S/N, and therefore the decision 
which frequencies to include, was computed from the com- 
bined data set. 

An interesting result is the large amplitude variability of 
two frequencies. The amplitude of 15.53 cd _1 varied from 
2.35 to 0.86 millimag from 2008 to 2009. The behavior 
of the new frequency at 31.60 cd _1 is puzzling. The fre- 
quency is essentially absent in 2008 and was present with a 
large amplitude of 2.1 millimag during 2009. Within each 
of the two observing seasons, the annual amplitude values 
did not change. Our attempts to identify the frequency as an 
instrumental artifact was unsuccessful: (i) Since the unre- 
duced raw data from 2009 also showed this frequency, it 
is not caused by our reduction scheme, (ii) The frequency 
is 2 cd _1 away from an orbital frequency (2f or b), but the 
observed amplitude is too high for an alias, (iii) Finally, 
the frequency value cannot be matched by a combination 
of another pulsation frequency with the orbital multiples. 
We conclude that the detection of the stellar 31.60 cd _1 fre- 
quency is probably real. 

We can now compare our frequencies with the known 
frequencies previously determined from ground-based pho- 
tometry. The agreement with the four frequencies deter- 
mined during the multisite Stephi campaign (Hernandez 
et al. 1998) is excellent. In fact, both studies find the same 
two frequencies to be dominant. The 34.23 cd _1 frequency 
found by Hernandez et al. (1998) (1.4 millimag amplitude) 
can be matched by 34.41 cd _1 (1.4 millimag amplitude) 
and a smaller-amplitude frequency at 34.20 cd _1 . Finally, 
the fourth frequency at 41.63 cd _1 agrees with our value 
of 41.65 cd _1 , although our amplitude of 0.24 millimag is 
much smaller. 

7 HD 73872 (= KW 375) 

The variability of this cluster star has not been known be- 
fore. A very brief variability test with photometric observa- 
tions covering only 2.8 hours showed constancy to 0.002 
mag (Breger 1970). However, the new data clearly show 
variability with small amplitudes. Table [4] lists the 18 fre- 
quencies detected by our MOST satellite data. 

8 Discussion 

The four Praesepe stars studied show a wide range of pul- 
sation behavior with frequencies from 0.76 cd _1 to 41.7 
cd _1 . The comparison of observed frequencies with pul- 
sation models (as well as individual stars with each other) 
requires the determination of the values of the pulsational 



constant, Q, for each frequency, i.e., the stellar mean den- 
sity needs to be known. The Q values are computed from 

logQ = -log/+0.51og^+0.1M bol +logT eff -6.456,(l) 

where the frequency, /, is measured in cd _1 . For S Set 
stars, narrowband uvbyj3 photometry is usually used to- 
gether with model-atmosphere calibrations to derive g and 
T e ff . Due to the uncertainty of the zero-points of the uvbyj3 
system and the presence of spectral lines, this procedure can 
lead to large uncertainties of up to 25% in the derived Q 
values. Fig. 7 of Breger (2000) demonstrates the problem. 
In this figure, a number of commonly used calibrations are 
used to derive a consistency mass, M, from the atmospheric 
parameters by 

M(xgL/T^. (2) 

A comparison of these consistency masses with the ex- 
pected evolutionary masses shows systematic deviations (up 
to 50%) as a function of temperature. Systematic deviations 
occur in all available model-atmosphere calibrations which 
we have tested. Keeping this warning in mind, we shall 
now calculate the Q values for the frequencies detected in 
the present study. The absolute magnitudes of the Praesepe 
stars are accurately known from the revised HIPP ARC OS 
parallaxes (Van Leeuwen 2009). We use the unpublished 
model-atmosphere uvby/3 calibrations by Villa (1998) and 
the senior author. This calibration uses c\ zero points to 
provide an agreement between the observed zero-age-main 
sequence (ZAMS) with the theoretical logg values along 
the ZAMS. This also minimizes the discrepancies between 
the computed consistency and evolutionary masses. We note 
that for A stars, our temperatures are close to those of the 
Moon & Dworetsky (1985) calibration, while the log g val- 
ues are, on the average, higher by 0.15. For the star BS Cnc, 
our value of 7600 K is in perfect agreement with the value 
published by Hernandez et al. (1998). 

Our derived values for T e ^ and log g are listed in Table 

m 

We can now normalize the observed frequencies to ac- 
count for the different mean stellar densities. We use Q = 
0.033d, which corresponds to the value for the radial fun- 
damental mode in S Set stars, as our reference value for 
Fig. [3] Furthermore, the scales on the frequency axes were 
chosen in such a way that all four panels have the same 
range in the pulsation constant, Q. This figure shows that 
the four stars exhibit different pulsation behavior, even after 
these frequency normalizations. Two stars (EP Cnc and BS 
Cnc) show high-amplitude peaks near the radial fundamen- 
tal mode, while the other two stars tend to pulsate in higher 
overtones. 

Several low frequencies outside the p-mode range were 
also detected, e.g., the 1.77 cd _1 peak for EP Cnc. Pre- 
vious earth-based telescope measurements have identified 
many low frequencies as combination frequencies, f^ — fj, 
where f^ and fj are two higher- frequency pressure modes. A 
beautiful example is the star 44 Tau (Breger & Lenz 2008). 
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Recent results from space missions such as Kepler have 
shown that low frequencies can also be a signature of hybrid 
pulsation, i.e., many S Set stars also pulsate with (7-Dor) 
gravity modes pulsation (Grigahcene at al. 2010). Further- 
more, stellar rotation commonly produces low-frequency 
peaks, both as f rot and 2f rot peaks (Balona 2011; Breger 
201 1). The four Praesepe stars have measured v sin i values 
in a narrow range, from 135 to 156 kms -1 . These values, 
together with the known absolute magnitudes and tempera- 
tures, make it possible to compute lower limits to the rota- 
tional frequencies. The lower limits range from 0.64 cd _1 
for BT Cnc to 1.6 cd" 1 for HD 73872. An interpretation of 
the observed low frequencies in terms of stellar rotation can, 
therefore, not be excluded at this point. 

Detailed asteroseismic interpretations are aided by 
mode identifications of a number of the detected frequen- 
cies. Such mode identifications can be obtained through reg- 
ular frequency patterns in the frequency spectra, multicolor 
photometry to give i values, as well as spectroscopic line- 
profile analyses. Regular frequency patterns cannot be seen 
in Fig. |3j the present data are single-color, and the ampli- 
tudes are too small for present-day line-profile analyses. 
Consequently, we are unable to provide mode identifica- 
tions at the present time. The present paper is only a first 
step of a larger MOST study of the Praesepe pulsators. We 
therefore defer the asteroseismic modeling to a later paper. 
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Table 3 Pulsation frequencies of BS Cnc for 2008 and 

2009 Table 4 Pulsation frequencies of HD 73872 for 2008 



Frequency Amplitude S/N Frequency Amplitude S/N 

2008 2009 ratio cd -1 millimag ratio 

cd _1 millimag 



fi 
f 2 
f 3 

u 

h 
fe 

h 
f 8 

fg 

fio 

fn 

fl2 
fl3 
fl4 
fl5 
fl6 
fl7 
fl8 



± u.uuz 

33.416 
35.981 
29.807 
33.542 
33.786 
26.086 
24.657 
30.685 
3.332 
32.215 
35.824 
31.016 
18.965 
27.241 
36.337 
31.947 
23.681 
31.714 



± u.uz 

3.03 
2.53 
1.52 
1.50 
0.88 
0.81 
0.76 
0.54 
0.47 
0.46 
0.46 
0.45 
0.39 
0.32 
0.27 
0.23 
0.22 
0.19 



29 
26 
21 
24 
17 
14 
15 
11 
8.6 
11 
10 
9.5 
9.6 
6.0 
6.1 
5.6 
5.7 
4.6 



35 
18 
23 
21 
18 
17 
15 
14 
14 
12 
10 
8.2 
8.0 
7.5 
7.4 
7.2 
7.3 
6.7 
6.3 
6.3 
6.2 
6.1 
5.9 
5.8 
5.8 
5.7 
5.3 
5.3 
5.2 
5.1 
5.1 
5.0 
4.7 
4.7 



' The formal errors for the listed frequencies are as small as 
0.0001 cd _1 . However, since annual alias values cannot be 
excluded, we list the value of the annual alias. 
2 ) The sum of these two frequencies is three times the orbital 
frequency, suggesting possible instrumental effects affecting the 
two peaks and their amplitudes. The agreement may be accidental, 
but some caution is advised. 



Table 5 Adopted atmospheric parameters for the four 
Praesepe stars 



Star 


T eff in K 


log# 




± 100 


±0.15 


BT Cnc 


7370 


3.45 


EP Cnc 


7950 


3.70 


BSCnc 


7600 


4.20 


HD 73872 


7750 


4.08 





± 0.0027 1} 


f. 


17 0363 


fo 


34 4088 




23.8118 


f 4 


15.5272 


f 5 


16.9347 


fe 


19.2807 


f? 


19.3379 


fs 


23.3644 


fg 


17.7465 


fio 


31.5988 


fii 


30.2671 


fl2 


6.38830 


fl3 


27.7331 


fl4 


36.6705 


fl5 


35.0716 


fl6 


34.1969 


fl7 


21.6585 


fl8 


26.6687 


fig 


27.5303 


f20 


37.0110 


f21 


36.8421 


f22 


21.1147 2) 


f23 


32.1866 


f24 


7.86910 


f 2 5 


21.4462 2) 


f26 


34.8510 


f27 


24.4267 


f28 


33.9294 


f2g 


35.1557 


f30 


41.6514 


f 3 l 


13.9000 


f32 


19.9478 


f33 


25.7046 


f34 


39.2526 



± 0.02 


± 0.04 


6.23 


5.72 


1.69 


1.15 




1 1 a 
1.13 


Z.J D 


U.oo 


1 1 Q 

L.Ly 


1 CM. 




n 87 

U.o / 


u.oo 


\j. i o 


U.OJ 


U.OJ 


U.oZ 




U.ZrU 


9 10 


n Q8 


U.jj 


n ^zl 

U. J^f 


U.JO 


u. / z 


u.UD 


^9 


99 


\J.DD 


U.J 1 




U.U/ 


0.33 


0.34 


0.38 


0.33 


0.38 


0.29 


0.27 


0.13 


0.25 


0.22 


0.26 


0.32 


0.40 


0.29 


0.27 


0.17 


0.24 


0.30 


0.23 


0.33 


0.28 


0.24 


0.39 


0.07 


0.22 


0.21 


0.22 


0.25 


0.47 


0.25 


0.23 


0.24 


0.17 


0.38 


0.21 


0.18 
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